FERMILAB-PUB-09-597-T 
MADPH-09-1550 



Model-Independent Searches for New Quarks at the LHC 

Anupama Atre^'^, Georges Azuelos^'^, Marcela Carena^'^, 
Tao Han^, Erkcan Ozcan^, Jose Santiago^ and Gokhan Unel^ 

^ Fermilab, Batavia U.S.A. 
^ Michigan State University, 
East Lansing U.S.A. 
^ University of Montreal, Montreal CA 

^ TRIUMF, Vancouver CA 
^ University of Chicago, Chicago U.S.A. 
^ University of Wisconsin, Madison U.S.A. 
University College London, U.K. 
^ University of Granada, Cranada Spain 
^ University of California, Irvine USA. 

Abstract 

New vector-like quarks can have sizable couplings to first generation quarks without conflicting 
with current experimental constraints. The coupling with valence quarks and unique kinematics 
make single production the optimal discovery process. We perform a model- independent analysis 
of the discovery reach at the Large Hadron Collider for new vector-like quarks considering single 
production and subsequent decays via electroweak interactions. An early LHC run with 7 TeV 
center of mass energy and 1 fb~^ of integrated luminosity can probe heavy quark masses up to 1 
TeV and can be competitive with the Tevatron reach of 10 fb~^. The LHC with 14 TeV center of 
mass energy and 100 fb~^ of integrated luminosity can probe heavy quark masses up to 3.7 TeV 
for order one couplings. 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) has finally started its quest to find the origin of 
electroweak symmetry breaking (EWSB). A natural realization of EWSB is likely to bring 
along new particles to cancel the ultraviolet sensitivity of the Higgs mass or to unitarize 
longitudinal gauge boson scattering in the absence of a light Higgs boson. A common 
occurrence among these new particles is new vector-like fermions that mix with Standard 
Model (SM) particles. Traditionally, vector-like quarks are considered to mix significantly 
only with the top sector, as the usual lore states that large mixings with lighter generations 
are excluded by electroweak precision or fiavor observables. This statement is however not 
necessarily correct, as cancellations among the effects of different types of new quarks can 
significantly alleviate the indirect constraints [iml. It was recently pointed out [4J that 
these cancellations can be the result of a symmetry |5| and naturally arise in certain models 
of extra dimensions or strong EWSB, which also allow for relatively light 3rd generation 
partners |6Nl3|. Motivated by this type of models, several studies of the LHC reach for 



top partners [14i-|l7l| and tau partners [18j have appeared recently in the literature. There 
are also some studies of new quarks mixing with first generation quarks. The discovery 
prospects of heavy down type vector-like quarks decaying to light quark generations in the 
context of the Eq GUT model I20I were studied in Refs. (21I . [22| for pair production mode 



and in Ref. [23|] for the single production mode. Most recently a more general analysis of 
the Tevatron reach for general new quarks with sizable couplings to the up or down quark 
was performed in Ref. j^. The goal of this article is to investigate the LHC potential in 
the search for general new heavy quarks exploring similar channels as in the previous study 
for the Tevatron. 

Heavy quarks can be produced in pairs via strong QCD interactions 

Qq, 99 QQ, (1) 

where Q generically denotes a new heavy quark. This is the most common and model- 
independent production mechanism for any states with QCD interactions. The heavy quarks 
can also be produced singly in fiavor-changing processes via electroweak interactions 

qq' ^ qiQ, (2) 

where V = W or Z gauge boson. Some of the current authors showed in Ref. [3] that the 
Tevatron reach to find new vector-like quarks with sizable mixings to first generation SM 
quarks is significantly better in the single production channel than in pair production. 

In this article we perform a theory-unbiased study of single production of new quarks 
with sizable mixing to the up or down quarks at the LHC. The large energy and luminosity 
of the LHC allow us to probe heavy new quark masses, even with rather moderate mixing to 
the first generation SM quarks, thus going beyond the need of exact cancellations to ensure 
compatibility with indirect constraints. We find quite encouraging results showing that an 
early run with 1 fb~^ of integrated luminosity at a center of mass energy of 7 TeV can be 
competitive with the Tevatron reach of 10 fb~^. The mass reach with a larger integrated 
luminosity of 100 fb~^ at a center of mass energy of 14 TeV can be up to 3.7 TeV for order 
one couplings. 

Amongst the many searches for new heavy quarks at collider experimentSjOne of the 



popular scenarios searched for is that of a sequential fourth generation quark j24j-|28|. Besides 
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the limits on their masses from direct searches at the Tevatron there are significant, recently 
updated bounds from electroweak (EW) precision data 25N27l|. The main difference between 
this model and the models under our consideration is that vector-like fermions, contrary to 
chiral fermions, decouple in the limit of large masses and therefore their effects can be 
made arbitrarily small, thus avoiding the severe constraints [295]. Phenomenologically, our 
theory parameterization and searching strategy are equally applicable to other heavy quark 
searches, as long as there is a sizable mixing with the light generations. 

The rest of the paper is organized as follows. In Sec. [Tll we present a general model 
parameterization for the heavy quarks and discuss the constraints on various models. In 
Sec. mil we calculate the heavy quark production cross sections at the LHC for the single 
production processes. In Sec. HVt we calculate the heavy quark decay via the charged and 
neutral currents and parameterize them by a few model independent parameters. In Sec. |Vl 
we explore the observability of the heavy quarks at the LHC via the single production 
mechanism with the energy options for 7 and 14 TeV. We present a discussion of our results 
in Sec. rvTland conclusions in Sec. IVIII 



II. GENERAL PARAMETERIZATION 



In this section we present a general parameterization for models with new heavy quarks. 
Let us consider new quarks labeled hj Q = X, U, D, Y with electric charges 
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(3) 



with masses rriQ and arbitrary couplings to the SM gauge bosons and first generation quarks 
parameterized by 



+ 



2cw ^ 



l^uU UrYUr + KdD dR'jf'DR 



+ h.c. . 



KuXUR-f^XR + KdudR-^^UR 



(4) 



This choice of electric charges exhausts the possibilities of new quarks mixing through gauge 
couplings with SM quarks. For simplicity we have only considered RH couplings, which is 
not a restriction since we are not using angular correlations in our analysis. Possible Yukawa 
couplings between U and u or D and d have not been written explicitly since they will not 
be used in our analysis as discovery channels. Finally, we have included only one new quark 
of each charge as multiple new quarks for each charge can be added trivially. 



A. Generic Constraints on the Couplings 

The couplings in Eq. (jl]) are subject to various constraints that can be classified as 
theoretical constraints and indirect experimental constraints. Theoretical constraints arise 
from the fact that HqQ originates from the mixing of SM quarks with new vector-like quarks 
and it is therefore bound by the unitarity of the matrices involved in such mixing. Generically 
we have 

K<1, 
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although this constraint could be stronger in particular models (see below). Indirect exper- 
imental constraints are much more model-dependent. They can arise at tree level through 
modification of the SM fermion couplings or at one loop, mainly through contributions to 
the oblique parameters 3^. The former crucially depend on whether cancellations among 



the contributions due to the mixing with different heavy quarks are present or not. In the 
absence of cancellations, the couplings HgQ can be quite constrained from low-energy exper- 
imental data. There is a generic correlation between the new couplings and corrections to 
the couplings of SM fermions, namely 

SgSM 2 /r-\ 
K , (5) 

gsM 

where qsm represents a generic SM coupling (see for instance Ref. [3l[). From EW precision 
observables, the SM couplings have been measured with a per mille precision. So, as a rule 
of thumb, we can expect a generic bound 

K < few X 10~^ (no cancellations). (6) 

As we will see below, cancellations among different contributions can however occur and 
make this bound completely disappear. One loop contributions to the oblique parameters, 
although also model dependent, are somewhat more robust. The resulting bounds on KgQ, 
which are typically milder than those in Eq. (jS]), can nevertheless be the main constraint in 
some models. 



B. Model Considerations 

As we have emphasized in the previous section, tree level experimental constraints on 
the values of KqQ are very model dependent and can be, for a given theoretical model, a 
lot milder than that of Eq. The reason is that there are models in which some of the 
KqQ can be sizable without modifying in an observable way the SM quark couplings (which 
are experimentally measured) and without inducing any observable flavor violation. This is 
the result of cancellations in the SM couplings due to the mixing of SM quarks with several 
different new quarks. These constraints or lack thereof can be better understood with a 
couple of illustrative examples as follows. 



1. Model I: Two degenerate vector-like quark doublets 

The flrst model, discussed in detail in Ref. [3], consists of two new vector-like quark 
electroweak doublets, with hypercharges 1/6 and 7/6 respectively, that mix only with the 
up quark, in the basis of diagonal Yukawa couplings for the charge 2/3 SM quarks. The 
particle content of these two doublets is two charge 2/3, one charge —1/3 and one charge 
5/3 quark. If the two doublets are degenerate, i.e. they have the same mass and Yukawa 
couplings to u before EWSB, their mixing with the up quark does not induce any observable 
correction to the SM quark couplings. This is due to the exact cancellation of two large 
contributions (arising from the mixing of the up quark with the charge 2/3 quarks in each 
doublet) with opposite signs. These large contributions are functions of KqQ and since the 
large contributions cancel, we can have large values of k^q without upsetting any SM quark 



4 



couplings. The only source of flavour violation in this case turns out to be suppressed by 
Cabbibo-Kobayashi-Maskawa (CKM) angles and by the ratio of the mass of the up quark and 
the mass of the new quarks, which provides enough suppression to keep flavor violation below 
experimental limits (details can be found in the Appendix of Ref. |^). At the loop level, 
the new quarks do however contribute to oblique parameters. In the degenerate case, the 
new sector is custodially invariant and the T parameter receives no correction. Interestingly 
enough, the contribution to the S parameter turns out to be essentially independent of the 
mass of the heavy quarks, for fixed values of Kuu- Note that for a fixed Yukawa coupling 
between the new quarks and the up quark, Kuu goes to zero in the large mass limit and 
therefore the vector-like quarks decouple in the heavy limit as they should. The dependence 
on the coupling Kuu is displayed in Fig. [H Strict constraints based on the contribution to 
the oblique parameters are difficult to impose without knowledge of other possible sectors of 
the theory that could also contribute to them. Nevertheless, assuming S < 0.2 (0.1) gives a 
pretty mild bound of < 0.75 (0.6) in the model under consideration. More details about 
the calculation of the S parameter can be found in Appendix \M 

0.4 
0.3 
^ 0.2 
0.1 
0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 1: One loop contribution to the S parameter, in the model with two degenerate doublets of 
hypercharge 1/6 and 7/6 and mixing only with the up quark, as a function of Kuu- The result is 
independent of the mass of the heavy quarks up to corrections 0{m'^/M'^). 

After EWSB one combination of the charge 2/3 heavy quarks couples to the up quark 
only through Yukawa couplings whereas the remaining three, of types U, D and X in our 
notation, have the following values of k, 

KuU = V2KuD = V2k,uX = Sr, (7) 

where sr is the sine of the corresponding mixing angle which depends on the particular 
values of the model parameters. Thus, we see that in this model, the generic bound from 
the unitarity of the mixing is saturated for Kuu < 1 and is a bit more stringent for KuD = 
i^ux < l/"\/2- We refer the reader to Ref. [4] for more details. 

Another similar model has two degenerate doublets of hyper charges 1/6 and —5/6, re- 
spectively, where the heavy quarks mix only with the down quark. In this case there would 
be one combination of charge —1/3 heavy quarks that couples to the down quark only 
through Yukawa couplings, while the other three quarks, U, D and Y in our notation, have 
large couplings to dn with the following constraints 

KdD = V2KdU = V2KdY = Sr, (8) 
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resulting in a theory upper bound of k^d < 1, t^du = i^dY < l/"\/2- As we have stressed 
above, these cancellations can be due to a symmetry and therefore be natural. 



2. Model II: One vector-like quark doublet 



In our second example no cancellations occur and SM couplings are modified due to 
the mixing with the heavy quarks. Stringent experimental constraints arise on the mixing 
and therefore on the values of KgQ (see Eqs. ([5]) — ([6])). The model consists of just one new 
electroweak doublet with hypercharge 1/6 that mixes only with the up quark. The particle 
content is one U and one D quark and the mixing with the up quark induces an anomalous 
ZurUr coupling, proportional to k^^, which is strongly constrained experimentally. A 



detailed analysis, using an updated version of the code in Refs. [32|, |33| gives a bound of 



^ 0.07. ^ The couplings to the heavy quarks can be computed as discussed in Ref. 



and Part 13 of Ref. 



(see Ref. [35[ for explicit expressions). Denoting the sine and cosine 



of the mixing angle by sr and cr, respectively, we have 

f^uU = SrCr, KuD = Sr. (9) 

The indirect bound on KuD, from the constraint on the ZurUr coupling, is therefore much 
more restrictive than the one coming from the unitarity of the mixing < 1- It also 
induces a very restrictive constraint on HuU ^ 0.07. Incidentally, there are also loop con- 
tributions in this model to the S and T parameters but the tree level constraints are much 
more restrictive. 



3. Model III: Nan- degenerate quarks 



An intermediate model between our two examples, with two new doublets which are 
not exactly degenerate allows for partial cancellations. For instance, using the code in 



Refs. 32|, |33j we find that Kuu ~ 0.25 is possible if it is within ~ 10% of \/2kux- Assuming 
similar values of the Yukawa couplings, this can translate to a 5 — 15% difference between 
the masses of different quarks, which might be experimentally distinguishable. 

Thus we can see that a large number of different possibilities can arise depending on the 
model but still sizable values of k^q cannot be excluded on a general basis. 



C. Generic Parameterization 



In order to consider all these possibilities, we do a model-independent study, imposing 
no a priori constraints on the values of KqQ. This way our results can be applied to any 
specific model, for which the relevant theoretical or experimental constraints will have to be 
taken into account. 



The code in Refs. [32, l33l assumes symmetry between the first two generations. The reported bound 
therefore apphes to a model with two new vector-hke muhiplets that mix equaUy with the first two 
generations. This can be considered a conservative but indicative bound in our model. 
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We can parameterize the coupling KqQ in a model-independent manner as 



l^qQ = {v/mQ)kqQ, 



(10) 



where the dimensionless parameter KqQ encodes the model-dependence and v ~ 174 GeV is 
the Higgs vacuum expectation value. This parameterization is useful because we generically 
have 



where A is some Yukawa coupling between SM and heavy quarks. Thus, in the limit that 
KqQ is relatively small so that the expansion can be well approximated by the first term, KqQ 
has a direct relation to a model parameter. It corresponds to a Yukawa coupling and it is 
naturally of order unity in a generic weakly coupled theory, and could therefore be up to 47r in 
strongly coupled theories. For larger values of KqQ non-linear effects become important and 
the dependence of kqQ on the model parameters becomes much more complicated. It should 
be noted that all the bounds we have obtained for the couplings, both the experimental 
ones and the ones coming from the unitarity of the mixing, apply to k and not k. Therefore 
heavier quarks can have a substantially larger value of k compatible with EW precision 
observables and unitarity of the mixing. 

III. HEAVY QUARK PRODUCTION 

In this article, we discuss the single production of heavy quarks via the process shown 
in Eq. ([2]). In order to keep our study widely applicable, we have factored out the model 
dependent coupling k"^, so that the cross section depends only on the corresponding Parton 
Distribution Functions (PDFs) and the kinematics (mass of the heavy quark produced). 
These model-independent cross sections will be denoted as cTq*" and aQ^ for single produc- 
tion of quark Q via charged current (CC) and neutral current (NC) interactions respectively. 
In the rest of the analysis, we adopt the CTEQ6L1 PDFs [36| for the parton distribution 
functions. 

The single production cross section as a function of the mass of the heavy quark is shown 
in Figs. [2] — m for various quark species at two different cm. energies, 7 and 14 TeV. 
We show separately the production cross section for each quark (solid black curve) and 
anti-quark (dashed black curve) species through charged (W*) or neutral (Z*) currents for 
single production. All the curves for single production correspond to a factorization and 
renormalization scale fir = fJ^R = mwimz) for CC (NC) processes. The scale dependence of 
our LO calculation is represented by colored bands that correspond to a scale variation up 
to fiF = fJ'R = riiQ. For comparison, the QCD cross section for pair production (independent 
of k) is also included in the figures and represented by dotted (black) curves. For pair 
production mode, we choose the central value of the factorization and renormalization scales 
tobe fip = = and the scale dependence of our LO calculation is represented by colored 
bands that correspond to a scale variation up to \/I/2 < jJ^F = I^r 1^ 2\/S. 

Fig. E] and Fig. [3] correspond to single production of D and U, respectively, with the 
left panels representing charged current channels and the right ones neutral current chan- 
nels. Although we have included all contributing channels for the production, the leading 




(11) 
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viQ (GeV) ruQ (GeV) 

FIG. 2: (a) Top left: cross sections (in fb) for singly producing the heavy quark D via CC 
interactions versus its mass mq at the LHC with y/s = 7 TeV in pp collisions, (b) top right: same 
as (a) but for NC interactions, (c) bottom left: same as (a) but for y/s = 14 TeV. (d) bottom right: 
same as (b) but for ^/s = 14 TeV. Solid (dashed) curves represent production of D{D) via the 
process pp — >■ jD{D). The cross section for pair production (which is independent of k) of heavy 
quark is shown as dotted (black) curve for comparison. The colored bands represent the variation 
in the leading order cross section due to the different scale choices as described in the text. 

contributions are from 

au — 7- jD, qa — )■ jU, 

ud ^ jU, qu ^ jU, (12) 

where q denotes a generic valence quark parton and j a light quark jet. Fig. HJon the other 
hand provides the production of X on the left and Y on the right. Due to their electric 
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FIG. 3: Same as Fig. [2] but for U. 

charge, they have only charged current channels and the leading contributions are from 

uu jX, dd jY. (13) 

Several interesting features regarding single production at the LHC can be observed from 
Figs. [2] — m To start with, let us recall the advantages of single production versus pair 
production as shown in Ref. jl]. The cross section for pair production falls off sharply due 
to phase space suppression and decreasing parton luminosity at large x values. Evidently the 
LHC with its substantially larger cm. energy will be able to probe higher quark masses in 
the pair production channel compared to the Tevatron. However, as at the Tevatron, single 
heavy quark production has the advantage of less phase-space suppression and longitudinal 
gauge boson enhancement of order m^/My at higher energies compared to pair production. 
Due to the participation of u, d valence quarks in the initial state with the coupling strength 
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niQ (GeV) mq (GeV) 



FIG. 4: (a) Top left: cross sections (in fb) for singly producing the heavy quark X via CC 
interactions versus its mass mq at the LHC with y/s = 7 TeV in pp collisions, (b) top right: same 
as (a) but for the heavy quark Y. (c) bottom left: same as (a) but for ^/s = 14 TeV. (d) bottom 
right: same as (b) but for ^/s = 14 TeV. The description of the curves is the same as Fig. [2] but 
for X and Y with only charged current interactions. 

given in Eq. (jlj), the cross section can be substantial and it falls more slowly for a higher 
mass. These effects can be easily observed from Figs. [2] — |H To summarize, the large cm. 
energy of the LHC coupled with the above advantages of single production make this an 
ideal process for discovery of heavy quarks up to very large mass values and very small 
couplings. 

Heavy quarks are produced at a much higher rate compared to anti-quarks with a larger 
difference for higher masses at the LHC. This difference in the rates of production of Q 
and Q is due to the difference in the PDFs of valence and sea quarks in the initial states 
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0.034 


0.0088 
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0.025 


0.0064 
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0.12 
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TABLE I: Cross sections in pb for EW single production of D, U, X and Y with k = 1 and 
mq = 900(1800) GeV at the LHC with ^ = 7(14) TeV. 



- heavy quarks couple to valence quarks whereas heavy anti-quarks couple to sea quarks. 
As we will see in later sections this feature will help us in constraining the electromagnetic 
charge of the heavy quarks. For jQ production valence quarks participate in the initial 
state and as the quark PDFs dominate at higher energies, the cross section falls slowly. Also 
note that, the relative sizes of the single production cross sections between different types 
of quarks are determined by the corresponding valence quark density in the initial state. 
To illustrate the features of single production explicitly, the cross sections for EW single 
production computed with MadGraph/MadEvent js^] are listed in Table |T] for all Q with 
rriQ = 900 (1800) GeV and k = 1 at the LHC with = 7 (14) TeV. These results reiterate 
the previous observations that (i) valence quarks dominate; (ii) u quark contributes more 
than d] and (iii) NC is weaker than CC. 

IV. HEAVY QUARK DECAY 

The singly produced heavy quarks will decay into jets and gauge or Higgs bosons through 
gauge and Yukawa interactions. In this article we are not considering the Higgs related 
channels so the relevant decay channels are 

X W+u, Y ^ W-d, (14) 

U -> W+d, Zu, D W-u, Zd, (15) 

and the charge conjugated ones. 

To perform a model-independent study, we classify processes according to their final 
states. In order to overcome the QCD background, we consider leptonic decays of the 
corresponding gauge boson. Therefore, the four final state configurations that could be 
experimentally distinguished from the SM backgrounds are: 

jj{w-^ru), jj{w+^i+u), jj(z^i+r), jjiz^uu). (le) 

Following this classification, we parameterize the corresponding cross section, under the 
narrow width approximation as 

cripp ^ nff) = S'^BviV ^ //'), (17) 
where V = W~^, W~ and Z with //' = i^u,l~v and plus vv, respectively. 
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All the model dependence is encoded in the 5*^ parameters which are nothing but the 
production cross sections times the corresponding branching ratios summed over all possible 
channels. These can be written as: 
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, ~2 CC 



(18) 



(19) 



(20) 



where cTq^'^^ are defined in Sec. [1111 and displayed in Figs. [2] — HI 

These parameters are what we are sensitive to from the experimental point of view. 
Their apparently complicated form is due to the high degree of model independence but 
can simplify greatly in specific models. For instance, assuming two degenerate doublets 
of hypercharges 1/6 and 7/6 that only mix with the up quark we have three quarks 
of charges 5/3, 2/3 and —1/3, all with the same mass mg, each one decaying with 100% 
branching ratio to W~^u, Zu and W~u, respectively. The fourth quark with charge 2/3, 
decays via the Hu channel and we do not consider that here. In that case we have 



S 



w- 



k^ n^^ 



uD'^ D + '^uX'^X 



CC 



S 



W+ 



k' 
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ni2 



with kuu = V^kuD = V^kux < ^q/v, where the inequality is the theoretical bound as 
discussed in Sec. Ill B 1[ In this specific case with two degenerate doublets it then implies 
the following upper bound on the corresponding 5*^, 



qw- ^ ^f^QVf^cc.cc 



CC\ 
X ) 



c:W+ ^ 1 ff^QY f^cc j_ ci 



CC\ 
D ) 



NC\ 
U > 



(22) 



The expressions for other models can be easily computed using the general expressions for 
the partial widths given in Appendix |Bl In the next section we describe the prospects for 
the observability of the new heavy quarks according to the signatures defined in this section. 



V. OBSERVABILITY OF THE HEAVY QUARK SIGNAL 

The classification of the possible final states that we have performed in the previous 
section allows us to perform a model-independent analysis of the LHC reach of new quarks 
which have sizable mixing with first generation quarks. We report on the results of such an 
analysis in the next few sections for single production. After a general discussion of the main 
kinematical features of the signal and background, we discuss the results of such an analysis 
for two different scenarios, one early run with 1 fb~^ integrated luminosity at y/s = 7 TeV 
and one longer run with 100 fb~^ of integrated luminosity at ^/s = 14 TeV. 
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FIG. 5: Illustration of the signal kinematics for a single heavy quark production. 



The final state channels of our current interest are 



(23) 



from Q decaying to a light-quark jet plus a W^^^ i^u) or ^(— ?■ uu), respectively. 

For the leptonic decay modes of the gauge bosons we consider i = e,^ and z/ = z/g, z/^, z/^. 
Although the inclusion of the r lepton in the final state could increase the signal statistics, for 
simplicity we ignore this experimentally more challenging channel. The relevant backgrounds 
to these processes are 

• W+ jets, Z+ jets with W, Z leptonic decays; 

• W~^W~ ,W^Z and ZZ with semi-leptonic decays; 

• single top production leading to W^b g; 

• top pair production with fully leptonic and semi-leptonic decays. 

An important aspect for our proposed search is the unique kinematics of single production 
and we discuss it first and present the detailed analysis next. A simple cartoon representation 
of the kinematic features is shown in Fig. O First, one of the two jets is from the heavy 
quark decay that makes it very energetic with a Jacobian peak in the transverse momentum 
spectrum near 



Using the px of the jets as a discriminant gives very good accuracy in identifying the jet 
coming from the decay of the heavy quark correctly, especially for high masses. Hence we 
identify the hardest jet (ji) in the event as the one from heavy quark decay. The second 
jet is associated with the t-channel exchange of a gauge boson. The coUinear behavior of 
the radiation implies a forward outgoing jet with its transverse momentum governed by the 




(24) 
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propagator's mass, typically around My/ 1 '2., and with a high rapidity. The pseudo-rapidity 
of the jet associated with the t-channel exchange of the gauge boson peaks at large values 
(typically \rj\ > 2). The soft nature of the radiated gauge boson makes the accompanying 
jet energetic. These features can be used for forward jet tagging to enhance the signal over 
the backgrounds. The forward tagged jet candidate is selected as the one with the highest 
deposited total energy and a minimum threshold in the pseudo-rapidity separation between 
the tag jet candidate and the decay products of the heavy quark, namely the reconstructed 
gauge boson and the highest pr jet, in addition to the basic cuts for jets during generation. 
These selection criteria accurately identify the forward tagged jet (jt) as the one associated 
with the t-channel exchange of a gauge boson. Note that there is also an s-channel diagram 
where the jet is not necessarily forward but the contribution of this channel is negligible. 

Similar to the hard jet, the W/Z from the heavy quark decay is also very energetic. 
This leads to rather coUimated final state leptons in the signal while those leptons in the 
background still tend to be back-to-back. The central hard jet and the energetic gauge 
boson W/Z are both from the heavy quark decay and are thus nearly back-to-back in the 
transverse plane. The separation (Ai?) of the central hard jet and the tagged jet is typically 
large for the signal as well, while that for the QCD background tends to present a coUinear 
singularity. 

Of most importance is the reconstruction of the mass peak for the resonant particles. 
For a heavy quark decay with Z — >■ in the final state, it is straightforward to form the 
invariant mass of the heavy quark mg = M(ji, Z) for the signal and backgrounds. For the 
final state with W ^ iu and Z — )■ uu, one can define a cluster transverse mass variable to 
be 



( \/ptw,z + Mw,z + PTn ) - {Ptw,z + PTn \ 



However, even for the W ^ C.v channel a full invariant mass can be constructed by including 
the W rest mass approximation (details are presented in Sec. IV Cp . Next, we present some 
details about the simulation. 

We have performed a detailed simulation of both the signal and all the relevant back- 
grounds using Madgraph/Madevent [s^] for event generation at the partonic level and 
PYTHIA jssi for parton showering with initial and final state radiation, as well as hadroniza- 
tion. Fast detector simulation was performed using Delphes [s^, version 1.8b, with the 
particle flow option turned off and the jets reconstructed using the anti-i^'T algorithm with 
a radius of 0.7 and a minimum transverse momentum of 20 GeV. All other parameters were 
set to default values that correspond to the ATLAS Detector. Missing energy was manually 
corrected for the muons in the event. 

The CTEQ6L1 PDFs j36| were used for both signal and background samples. For the 
signal we set the factorization and renormalization scales to be /zf = /i/? = mw{fnz) for 
CC (NC) processes. Due to the nature of the vector boson fusion, it is known that the 
natural choice of the factorization scale is around the mass of the propagating particle 



[4d . |41| . Varying the scale within some reasonable range could result in an uncertainty 
of about 10%. For the background processes the renormalization and factorization scales 
are set to jJip = I^r = = '^PrU))- To avoid double counting of the phase space 

overlap between the parton shower and matrix element parton generation, the modified MLM 



matching procedure |42| with the showerkt method [43| was applied for the background 
processes with multijets in the final state. The value of the parameter xqcut, which sets the 
minimum Kt measure of partons in the matrix element generation, was set to 30 GeV. We 
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TABLE II: Cross section (in pb) for the various background processes with the acceptance cuts in 
Eq. (p5]) . A requirement that priJi) > 40GeV was applied in ah cases, except where indicated. 



Background 


y/s =7TcV 


yfs =14TeV 


(W ^ tv) + 1,2,3 jets 


1070 


2680 


{W -^iv) +2 jets 






VTi^^ ) > 400 GeV 




2.76 


{Z vv) + 1,2,3 jets 


402 


1206 


{Z It) + 1,2,3 jets 


97.8 


247 


w± pi/T J- 1 2 iets 


7.83 


19.0 


Zii + 0,1 jets 


1.64 


4.18 


Wf^ Zu + 0,1 jets 


0.429 


1.05 


Zu Zjj +0,1,2 jets 


0.291 


0.689 


single top (t — )• £^vb) + 1,2 jets 


7.24 


28.8 


tt fully leptonic 


3.73 


21.1 


tt semi leptonic 


20.4 


121 



summarize below the minimum cuts imposed in reconstructing the physics objects for our 
analysis. 

Pt(ji) > 40 GeV, p'^ > 20 GeV, pi > 15 GeV, 

\rij\ < 5, \rie\ < 2.5, AR{j£) > 0.4, AR{ii) > 0.2. (25) 

Various background cross sections are given in Table [ITl Finally QCD di-jet samples 
at various hard-scattering momenta {p±) were generated with Pythia to be used as back- 
grounds to the Zi,i^ + 2 jets analysis. The cross-section for the 280 GeV< p± <560 GeV range, 
which is found to be the most relevant in terms of mimicking the signal, is 1.82 nb at 7 TeV 
collision energy. 

It is evident from Table [ITl that W/Z +iets is the leading background. However there is 
also a sizable contribution to the total background from ti production. While the cuts we 
describe later are sufficient to extract the signal from the large backgrounds, we list below 
some further improvements that can be made in order to reduce the large tt background. For 
the case of CC decay modes of the signal, we can impose a veto on a second isolated lepton 
to reduce the ti background from the fully leptonic decay mode. For the NC decay mode of 
the signal with leptons in the final state, we can veto events with large missing energy. To 
further reduce tt events, we can veto events with 6-tags (this has virtually no effect on our 
signal). The semi-leptonic decay mode can be reduced further by vetoing events where any 
two jets reconstruct a VT-boson or if any three jets reconstruct the top quark. These simple 
vetoes, together with the rest of the cuts allow us to reduce the ti background to negligible 
levels. In the next three subsections we present the CC and NC analyses making use of the 
unique kinematic features discussed above. 
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For the rest of the analyses, we consider a variety of heavy quark masses and study the 
reach at two LHC energies, as given below 



mg = 300, 600, 900 GeV, = 7 TeV; 

mg = 900, 1800, 2400 (or 3000) GeV, = 14 TeV. (26) 

A. Neutral Current Channel with Z-Decays to Charged Leptons 

In this section we study the case where the heavy quark decays into a jet and a Z 
boson, with the Z boson decaying leptonically into electrons and muons [l^ = e^, /z^) and 
investigate the prospects of the LHC in this channel. Since the Z'-boson's invariant mass can 
most easily be reconstructed using the charged leptons, the analysis starts with the selection 
of those. The basic acceptance cuts are listed in Eq. fl25|) . The invariant mass reconstructed 
using the two leptons is required to be within 25 GeV of the Z-boson mass peak to remove 
any possible fakes and misidentifications. Additionally, to select the boosted Z bosons from 
the signal, the azimuthal angle between the two leptons is required to be smaller than a 
maximum, the value of which is optimized for the best signal significance. 

The event is also required to contain at least two jets isolated from leptons and other jets, 
the criteria for which are listed in Eq. fl2^ . The jet with the highest is selected as the 
one coming from the heavy quark decay. The tag jet candidate is selected as the jet with 
the highest deposited energy and a minimum threshold on the pseudorapidity separation 
between the tag jet candidate and the decay products of the heavy quark, namely, the 
reconstructed Z boson and the jet with the highest pt- To take advantage of the forward 
nature of the tagged jet, thresholds optimized for best signal significance are applied for the 
pseudorapidity of the tagged jet. The efficiencies of these selection criteria in the order they 
are applied for various signal and background samples are shown in Table IIIII for the 7 TeV 
and 14 TeV cm. energy scenarios. 

The final heavy quark (Q) candidate's invariant mass is formed using the reconstructed 
Z boson candidate and the jet with the highest transverse momentum. The results of 
this are shown in Fig. [6t^a) for a heavy V type quark with mass of 600 GeV for ^/s = 
7 TeV. The black, blue and red histograms represent background, signal and signal plus 
background respectively. In order to extract the signal significance, a series of successive 
fits are performed to the signal+background sample. The procedure starts with a Landau 
functional fit which should describe the overall behaviour of the data set. The values of the 
parameters obtained from the Landau fit are used as initial values in a second fit which uses 



a Crystal Ball function |4J] . In the absence of signal the Crystal Ball function is expected 
to describe the sample behavior in an ideal way. The obtained curve is shown with a dashed 
black line in Fig. |6]^a). If the signal hypothesis is true, some mismatch between the fitted 
function and the histograms is expected. Therefore a new fit is performed to investigate the 
signal hypothesis. The fit function is the sum of two components: a Crystal Ball term to 
represent the background and a Breit-Wigner or a Gauss term term to represent the signal 
resonance. The selection between Breit-Wigner and Gauss terms is made based on the mass 
of the sought heavy quark: up to 1 TeV, the signal is narrow, thus best represented by the 
former and after 1 TeV, the signal distribution becomes fatter thus compatible with the 
latter. The obtained curve is shown with a red solid line in the invariant mass distributions 
of Fig. El^a). The background and signal components of this function are also shown on 
the same plot as dashed red and solid magenta lines, respectively. Therefore a good match 
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TABLE III: Selection criteria and the efficiencies (in percent) for the signal (emg ) and background 
(ebg) for the case of the neutral current channel with the Z boson decaying to electrons or muons. 
Note that in evaluating esc all the relevant backgrounds have been included. 





LHC 7 TeV 


LHC 14 TeV 


selection 


limit 


£600 




limit 


cgoo 




limit 


£2400 




min ^ of e//i 


2 


92 


78 


2 


95 


78 


2 


98 


54 


\Mir-Mz\ < 


25 


95 


87 


25 


95 


87 


25 


54 


40 


|A0f/| < 


0.9 


78 


26 


0.7 


86 


18 


0.3 


82 


6 


min ^ of jets 


2 


92 


67 


2 


92 


66 


2 


93 


90 


|A?7jijJ, \^Vz,n \ > 


1.4 


80 


40 


1.2 


86 


47 


1.5 


87 


48 


Pt > 


180 


89 


37 


250 


91 


16 


610 


84 


13 


En > 


180 


83 


54 


220 


83 


46 


380 


86 


51 


Vh > 


1.3 


97 


96 


1.8 


90 


87 


0.8 


100 


99 


overall 




36 


0.92 




41 


0.25 




26 


0.034 




FIG. 6: Reconstructed heavy quark distributions at the LHC with y/s = 7 TeV. (a) left: invariant 
mass of the heavy quark of U type and mass 600 GeV decaying leptonically via NC interactions. 
See text for a description of all the curves; (b) right: transverse mass of the heavy quark candidates 
of D type with mass 350 and 400 GeV decaying invisibly into neutrinos via NC interactions (after 
all the selection criteria have been applied). 

between the signal (solid blue histogram) and its estimation (solid magenta line) shows the 
power of the method described above. 

Finally, we mention some details about the signal simulation sample. We have simulated 
signal samples for a D type quark for the 7 TeV analysis and a U type quark for the 14 TeV 
analysis with k = 1 and Br{Q — )■ qZ) = 100%. While the efficiencies we list in Table UTTl are 
from one quark sample for each cm. energy, we have checked that there is no significant 
change in the efficiencies while using the other quark sample at the same energy. Note that 
the choice of the quark sample and the specific model parameters used does not affect the 
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model-independent quantity inferred from an excess of signal events over the background 
at a given luminosity. The choice of the quark sample does however affect the luminosity 
required for evidence or discovery of such a quark and the constraints on the various model 
parameters inferred from it. 



B. Neutral Current Channel with Z-Decays to Neutrinos 

In this section we study the case where the heavy quark decays into a jet and a Z 
boson, with the Z boson decaying invisibly into neutrinos. An accurate reconstruction of 
the neutral channel in the Z,^,^+jets final state is beyond the scope of this paper, largely 
due to difficulties in simulating the measurements of large missing transverse energy. Even 
with detailed GEANT4 simulations of the detectors, fake E^^^^ from event mismeasurements 
(due to non-Gaussian tails in the detector jet responses) is unlikely to be determined with 
confidence before studies are done with control samples in actual LHC data. However, this 
final state can still offer some supplementary evidence in the heavy quark searches as we show 
in this subsection. As an example, we consider the discovery potential for two relatively-low 
quark masses (350 and 400 GeV) with 1 fb~^ of data at 7 TeV cm. energy and describe the 
analysis in detail below. 

After all the selection criteria (as listed in Table HVT) . the dominant background is found 
to be the irreducible background from the Standard Model Zj^^-l-jets processes. A subdom- 
inant contribution is from the PV^i^-l-jet events, where the charged lepton falls outside of the 
detector acceptance or is lost in the reconstruction. We also consider semileptonic ti and 
jet-associated single top as other sources of background and find them essentially negligible. 

As can be seen in the table, the requirements on |A0| between the two leading jets and 
between any of the three leading jets and the missing transverse momentum appear to play 
no significant role in separating the signal from the main sources of background. However, we 



prefer to apply these cuts, since in Ref. [45j they have been shown to suppress backgrounds 
from QCD processes. All the cuts are tested on Pythia di-jet samples generated at various 
hard scattering momenta, and while we fail to reproduce the results from this reference due 
to limitations of the fast simulation, we observe that the cuts on the tag jet also provide 
an extra handle in reducing the QCD backgrounds to negligible levels. Note that for this 
channel, the tag jet also turns out to be the jet with the second highest transverse momentum 
i.e. it = 32 and the tag jet is denoted as j2 in this section. 

Attributing the total missing momentum in the event to the Z boson, the transverse mass 
of the heavy quark candidate is reconstructed as shown in Fig. [6t^b). While it is clear that 
the signal and background distributions have distinct shapes and therefore the signal could 
be extracted using a fitting technique, we conservatively choose to do a simple counting 
measurement for this final state. Since most of the background is from the Zy,^+]ets and 
to a lesser extent from the ly^j^-l-jet events, there are natural control samples obtained by 
reverting the charged lepton veto done at the beginning of the event selection. For example 
-^f£+jets events can be used to determine the amount of Z,j„+]et?> background. Therefore for 
an estimate of the expected significance, we calculate s/Vl/, where s is the expected number 
of signal events and b' includes the expected uncertainty in the background estimation. For 
simplicity, if we assume b' = b x BR{Z — > vv) / BR{Z ££), where i = e, fi and b is the 
expected total number of background events, we find that the significance for the 350 (400) 
GeV quark signal is about 6.0 (4.2)(j with 1 fb~^ of data. 

Next we give some details about the signal simulation sample. We have simulated signal 
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TABLE IV: Selection criteria and the efficiencies (in percent) for the signal and background for 
the case of the neutral current channel with the Z boson decaying to neutrinos. Note that j2 here 
is the same as jf 





LHC 7 TeV 


selection 


limit 


£350 


£400 






# of e//i < 


1 


99.5 


99.4 


99.5 


10 


Pr > 


100 


73 


81 


13 


11 


min # of jets 


2 


92 


92 


62 


65 


> 


100 


91 


94 


63 


59 


P'^ > 


40 


77 


80 


58 


52 


\Vj2 1 > 


1.0 


74 


73 


51 


57 




3.0 


45 


46 


15 


19 


E,. > 


100 


99 


99 


97 


99 




2.7 


86 


84 


82 


87 




0.5 


79 


78 


73 


59 


overall 




10.5 


12.1 


0.13 


0.01 



samples for D type quarks with k = 1 and Br{D — )■ dZ) = 100%. From the studies of 
the neutral current channel decaying to charged leptons and the charged current channel 
decaying leptonically, we have found that there is no significant change in the efficiencies 
while using a U quark sample. Hence, using the same efficiencies as listed in Table HVl for a 
U type quark, which has a larger production cross section compared to the D type quark, 
would have resulted in a smaller luminosity required for evidence or discovery. For simplicity, 
we only give the results from the D type quark as an illustration that this channel is a viable 
mode for evidence or discovery, without loss of generality, as discussed in the previous section. 

C. Charged Current Channel 

In this section we study the case where the heavy quark decays into a jet and a W boson, 
with the W boson decaying leptonically W ^ iu {£ = e, and investigate the prospects 
of the LHC in this channel. Once again, we adopt the mass parameters as in Eq. (126|) . and 
describe the analysis in detail below. 

The simplest way of reconstructing the W boson candidate in the charged channel is via 
its leptonic decays and hence the reconstruction algorithm starts with the selection of events 
containing at least one electron or muon {£) with the selection criteria described in Eq. 
In the case where both flavours of leptons are available, the lepton with the higher transverse 
momentum is selected as the W decay product candidate. After selecting events with a 
minimum missing transverse energy {^), the four momentum of the W boson is obtained 
using the information from and and including the W rest mass approximation. The 
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two-fold ambiguity in the solution of the W boson longitudinal momentum is handled by 
choosing the most central (minimum \r]w\) solution. We find that, while this method is only 
slightly better than a random choice at very high quark masses, it reconstructs the true 
momentum about 65% of the time for low heavy quark masses, and improves the resolution 
of the reconstructed heavy quark mass. The azimuthal angle between the charged and 
neutral lepton is required to be smaller than a threshold (again optimized for each mass 
point for the best signal significance) to select the boosted W bosons. 

The event is required to contain at least two jets isolated from leptons and other jets, the 
criteria for which are listed in Eq. f l25|) . The jet with the highest transverse momentum is 
assumed to be the decay product of the heavy quark. We identify the tag jet as the jet with 
the highest energy deposition in the event that is at least a minimum |A?7| apart from both 
objects that we use to reconstruct the heavy quark candidate i.e. the tag jet has to be well 
separated from both the reconstructed W boson and the highest pt jet. 

We introduce a new variable called "sum px" which is defined as the absolute value of 
the vector sum of the tagged jet and the reconstructed heavy quark transverse momenta to 
enhance the signal over backgrounds, 

(27) 

At the leading order, this variable is zero for the signal but detector effects and initial 
state radiation smear this distribution. However the distribution peaks at very small values 
for the signal while it is typically large for busy background events, in particular top pair 
production. We require the to be smaller than a threshold to enhance the signal over 
the backgrounds. This cut reduces the top pair background significantly for the 14 TeV 
scenario while also having a reasonable impact for the 7 TeV analysis. The efficiencies 
of these selection criteria in the order they are applied for various signal and background 
samples are shown in Table IVl for the 7 TeV and 14 TeV cm. energy scenarios and different 
mass values of the heavy quark. The superscript in the efficiencies of the backgrounds (for 
e.g. i in e^g.) correspond to ly+jets, diboson, top pair and single top for i = 0,1,2,3 
respectively. 

The leading jet and the W boson are combined to obtain the reconstructed invariant 
mass of the heavy quark candidate. The procedure to extract the signal significance is the 
same as described earlier in Sec. IV A[ The invariant mass distribution for a heavy quark 
of D type and mass 900 GeV decaying leptonically via CC channel is shown in Fig. [71 
In Fig. [Tl^a) we show the reconstructed heavy quark where the charge of the lepton is not 
identified and in Fig. El^b) when the charge of the negatively charged lepton coming from the 
decay of the heavy quark is identified. When we do not identify the charge of the lepton, we 
have to include signals from both quark and anti-quark and the background would include 
the contribution from both ly^+jets. This method adds a small component to the signal 
cross section as the anti-quark cross sections are much smaller than quark cross sections 
while adding significantly to the backgrounds as positively charged W bosons are produced 
preferentially at a pp collider such as the LHC This reduces the signal significance. However 
identifying the charge of the negative lepton enhances the signal over the backgrounds and 
reduces the luminosity required for discovery by almost half. For the heavy quarks decaying 
into positively charged leptons (say U and X), again the anti-quark contributions add very 
little to the signal but the reduction in the background is much smaller as are produced 
preferentially at the LHC. However, there is still some improvement in the signal significance 
by identifying the charge of the positive lepton though not as much as in the case of the 
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FIG. 7: Reconstructed heavy quark distributions for a quark of type D and mass 900 GeV decaying 
via CC interactions at the LHC with ^/s = 7 TeV. (a) left: no charge identification of the charged 
lepton has been made; (b) right: with charge identification of the negatively charged lepton. See 
Sec. IV Al for a description of all the curves. 

negative leptons. The effect of identifying the charge of the lepton (for both positive and 
negatively charged leptons) is shown in Table |Vl 

We show some sample distributions for a few of the relevant kinematic quantities which 
are unique to single production of heavy quarks and enable us to separate the signal from 
the background in Fig. |H1 The pr distribution of the hard jet and the energy distribution of 
the forward tagged jet are shown in Fig. IHI^a) and Fig. MJo) respectively. The difference in 
pseudo-rapidity between the hard jet and the forward tagged jet is shown in Fig. ^c) and 
the azimuthal angle separation between the leptons {p^ and '^t) is shown in Fig. |8]^d). All 
distributions in Fig. [8] are normalized to 1 pb~^ integrated luminosity. 

We mention some details about the signal simulation sample. We simulate signal samples 
for both D and U type quarks with k = 1 and Br{Q — )■ qW) = 100%. The efficiencies we list 
for the various cuts (except the charge identification cut) in Table |V] are from the D quark 
sample. However, we have checked that there is no significant change in the efficiencies 
while using a U quark sample, except for the cut on Ej^ where a U type quark has a higher 
signal efficiency compared to the D type quark. The reason is as follows: from Eq. f ll2p we 
see that the dominant contribution comes from the case where the tag jet is a SM d quark 
for D type quark and a SM u quark for the U type quark. Due to the PDFs, the u quark 
has a slightly harder energy distribution compared to the d quark, which results in a higher 
efficiency for the tag jet in the case of the signal with U type quark compared to the D type 
quark. While the numbers we list in Table |V] are for the D type quark we use the actual 
values for the U type quark when extracting the signal and background events for evidence 
or discovery reach. 

In the last few rows we show the efficiency when the charge of the lepton coming from the 
decay of the heavy quark is identified. For Qi = —1 we use the D type quark sample and for 
Qe = +1 we use the U type quark sample. Note that the efficiency for identifying a lepton 
correctly as positively or negatively charged is the same irrespective of origin of the lepton. 
In fact in our analysis the efficiency of identifying the charge of a lepton itself is essentially 
100% (see Sec. |Vl]for details). However from Table IVl it is evident that a U type quark has 
a slightly higher efficiency in identifying the charge of the lepton. The reason is that we 
start with a Q + Q sample and then identify the charge of the positive or negatively charged 
leptons in this sample. Given that U type quark has a smaller anti-quark component (which 
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TABLE V: Selection efficiencies (in percent) for signal (emg) and background (e^g) samples for the 
case of the charged current channel with leptonic decays into electrons and muons. The superscript 
in the efficiencies of the backgrounds (for e.g. i in €^^q) correspond to VF+jets, diboson, top pair 
and single top for i = 0, 1, 2, 3 respectively. 





LHC 7 TeV 


LHC 14 TeV 


selection 


limit 


£600 






^BG 




limit 


£900 


^BG 


A 


^BG 


''EG 


limit 


£3000 


£bg 


''EG 


'^EG 


£bg 


min ^ 
of e//i 


1 


96 


91 


91 


94 


92 


1 


97 


91 


91 


94 


92 


1 


99 


98 


90 


96 


96 


fjT > 


30 


94 


76 


70 


80 


80 


50 


92 


57 


35 


54 


52 


110 


96 


62 


8 


39 


65 


V solution 




82 


80 


74 


73 


74 




81 


80 


70 


71 


71 




67 


78 


76 


77 


77 


|A0<!,^| < 


0.8 


86 


41 


13 


31 


19 


0.6 


91 


29 


13 


25 


14 


0.5 


99 


84 


50 


63 


71 


min ^ 
of jets 


2 


91 


63 


89 


100 


95 


2 


90 


63 


90 


100 


95 


2 


90 


96 


95 


99 


95 


\^'n3un \ > 


1.8 


68 


28 


20 


30 


28 


2 


65 


24 


17 


25 


23 


2.3 


72 


31 


25 


32 


30 


Pt > 


200 


82 


25 


26 


18 


15 


250 


90 


17 


24 


14 


11 


450 


98 


67 


53 


24 


70 


En > 


200 


78 


61 


60 


71 


62 


300 


74 


52 


47 


60 


52 


270 


95 


72 


73 


79 


70 


\Vn 1 > 


0.6 


100 


100 


100 


100 


100 


0.6 


100 


100 


100 


100 


100 


0.6 


100 


100 


100 


100 


100 


Pt 


100 


94 


76 


56 


47 


47 


100 


93 


74 


56 


41 


40 


100 


74 


48 


38 


18 


18 


Qe = -1 




85 


31 


39 


49 


28 




88 


31 


50 


47 


36 




93 


29 


45 


45 


29 


overall 




20 


0.15 


0.037 


0.15 


0.031 




21 


0.036 


0.014 


0.035 


0.0089 




26 


0.82 


0.041 


0.088 


0.25 


Qe = +1 




90 


69 


61 


51 


72 




94 


69 


50 


53 


64 




97 


71 


55 


55 


71 


overall 




24 


0.32 


0.059 


0.16 


0.079 




25 


0.083 


0.014 


0.039 


0.016 




29 


1.99 


0.050 


0.11 


0.60 



gives the wrong sign leptons) of the total cross section compared to the D type quark, we 
get a shghtly higher efficiency for charge identification for the U type quark. 

Finally, note that the choice of the quark sample and the specific model parameters 
used does not affect the model-independent quantity inferred from an excess of signal 
events over the background at a given luminosity. The model-independent parameterization 
discussed in Sec. IIVI was designed such that the effects of k and the branching ratio of the 
heavy quark decay mode could be factored out. The choice of the quark sample does however 
affect the luminosity required for evidence or discovery of such a quark and the constraints 
on the various model parameters inferred from it. 
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FIG. 8: (a) Top left: pT of the hardest jet. This jet is used in reconstructing the heavy quark mass; 
(b) top right: total energy of the forward, tagged jet (c) bottom left: difference in pseudo-rapidity 
between the hardest jet and the tagged jet (d) bottom right: azimuthal angle between pt{^) and 

D. Identification of tiie Heavy Quark 



In this section we discuss the interesting possibility of determining the identity of the 
heavy quark produced. In pp collisions such as at the Tevatron, this is rather straight forward 
as demonstrated in Ref. Defining the proton beam as the +z direction, the forward 
(backward) tagged jet should be correlated with the heavy anti-quark (quark) production. 
Moreover, this can be used as an indication for down-type or up-type heavy quark production 
by identifying the electromagnetic charge of the lepton from the charged current decay. For 
e.g., an event with a backward tagged jet and a positive (negative) lepton would indicate 
production of U (D) heavy quark. Similarly, an event with a forward tagged jet and a 
positive (negative) lepton would indicate production of D (U) heavy quark. However, note 
that one cannot distinguish between a charge 2/3 and charge 5/3 quark (or a charge —1/3 
and charge —4/3 quark) by this procedure. 

At a pp machine such as the LHC, the situation is somewhat different. First of all, 
as seen from the total cross section curves in Figs. ([2]) - (j3]), a heavy quark production is 
substantially larger than that for an anti-quark, due to the u, d valence quark dominance. 
There is a further subtle difference between the production of a heavy quark and an anti- 
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quark. A valence quark parton has a harder distribution in the energy fraction x, while an 
anti-quark parton distribution is a lot softer. Consequently, a heavy quark produced will 
be more strongly boosted. Finally, the leptonic decay of a heavy quark via charged-current 
may still be utilized to specify the electric charge of U oi D. 

Another way to exploit the difference in the cross section of a heavy quark and anti-quark 
to gain sensitivity to the charge of the heavy quark produced is described below. We define 
a charge asymmetry which will be independent of the couplings k or of the branching ratios: 

A = ^^±^, (28) 

where N± is the cross section measured for the W^j channel. For example, at ^/s = 7 TeV, 
based on the values given in Table [H one finds A = (—0.77,0.84,0.95, —0.41) for quarks of 
mass 900 GeV and of type (D,U,X,Y). 

Note however, that both of the methods described above work best when the quarks are 
non-degenerate and the resonances can be distinguished experimentally. In that case, the 
charge asymmetry and/or the rapidity of the reconstructed quark along with the charge of 
the observed lepton could serve to distinguish between the quarks with different electromag- 
netic charges. When the states are exactly degenerate (equal mass and equal couplings) 
these methods have a contamination from the anti-quarks from a few percent up to few tens 
of percent depending on the model and the heavy quark under consideration. For example, 
in the case of doublets with hypercharges 1/6 and 7/6, upon observing a 5a excess in the 
channel for a mass of 900 GeV, one would have a contamination from the X quark 
~ 5% in the case of degenerate bidoublets. If the quarks were non-degenerate with cou- 
plings also varying by 10%, then the contamination is (9(10%). If the observation were to 
be in the channel for the hypercharge 1/6 and -5/6 doublets, the contribution from 

the anti-quarks increases to ~ 20% for the degenerate case and further goes up to (9(30%) 
for close by resonances with slightly unequal couplings. In summary, the methods described 
for charge identification are best applied for the case of non-degenerate quarks and caution 
must be exercised when using the same methods for degenerate (and very close in mass and 
coupling) quarks. 

VI. RESULTS 

In this section we will consider the current constraints from direct searches, the sensitivity 
of the search for new heavy vector-like quarks at the LHC, the constraints on some example 
models at the LHC and finally discuss some systematic effects. 



A. Current Constraints from Direct Searches 

The current bound from direct searches at the Tevatron experiments is 

niQ > 335 (268) GeV at 95% C.L. (29) 
for heavy up (down) type quarks produced in pairs and decaying via CC (NC) interactions 



I46l . 14 7j . Recently both CDF and DO experiments at the Tevatron did a search for heavy 



quarks with sizable couplings to first generation SM quarks in the single production channel 
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4i,|49|. They presented their results in a model-independent fashion in the plane of mass 
of the heavy quark vs cross section times branching ratio (a parameter similar to S^). From 
these analyses the heavy quarks have limits as given below. 



CDF : a{pp ^ qQ) x Br{Q ^ qW) < 240 fb for mg = 530 GeV 
DO : a{pp qQ) x Br{Q qW) < 40 fb for tuq = 690 GeV 
DO : a{pp qQ) x fir(Q ^ qZ) < 120 fb for tuq = 550 GeV, (30) 

where Q is a generic heavy quark. The charge of the lepton was not identified in these 
analyses and hence the correlation with our parameters is slightly different and is given 
by 

a{pp qQ) xBr{Q^ qW) = S^' + S^^ for CC, 
a{pp qQ) X Br{Q ^ qZ) = for NC. 

The Tevatron being a pp machine, heavy quarks and anti-quarks are produced in equal 
numbers with the result that = . The constraints from the Tevatron (on + 
) can hence be translated into constraints on S^"^ or . 
While it is very useful to have model-independent limits due to the wide applicability of 
results, it is also illustrative to interpret the results in the case of a particular model. For 
example, interpreting these results in the context of degenerate doublets (with hypercharges 
1/6 and 7/6) with kuo = 1 and decaying 100% via CC currents, the results from CDF exclude 
heavy D type quarks up to mass 530 GeV while the DO constraints are more stringent and 
exclude heavy D type quarks up to 690 GeV. In the same model for a U type quark with 
kuu = and decaying 100% via NC currents, the sensitivity is to a 550 GeV quark. Note 
that in the context of this model, since all the quarks are degenerate, a U type quark of 
mass 690 GeV decaying via NC is also excluded. However in other models the interpretation 
can be quite different. 



B. Sensitivity at the LHC 

We have simulated signal and background events for two different scenarios: an early run 
with ^/s = 7 TeV and an integrated luminosity of 1 fb~^, and a longer run with ^/s = 14 
TeV and 100 fb^^ integrated luminosity. The estimated number of signal and background 
events are extracted using the fitted function to the total number of events as described in 
Sec. IV A[ The statistical significance is calculated as [soj] 



S=^2x[{s + b)\n{l + ^)-s], (31) 

where s{b) is the number of signal (background) events determined from the Breit-Wigner 
(Crystal Ball) term of the fitted function by integrating over 2a of the most probable value 
of the Breit-Wigner component. The estimated number of signal and background events 
extracted from the fits to the total number of events for the LHC with ^/s = 7 TeV and 14 
TeV are presented in Table. I VI I and Table. IVIII respectively. For each case, the number of 
signal and background events giving rise to 3a evidence and 5a discovery (in parentheses) are 
also given. In the same tables are also given the luminosities required for 3a evidence and 5a 
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TABLE VI: The estimated number of signal and background events extracted from the fits to the 
total number of events for the LHC with y/s = 7 TeV. For each case, the values for 3a (and 5a in 
parentheses) statistical significance is also given. 



Channel : £+£- 2j 


mg (GeV) 


300 


600 


900 


Signal events 


10.7 (29.2) 


7.10 (20.3) 


6.80 (18.9) 


Background events 


9.30 (25.4) 


3.70 (11.0) 


3.30 (9.20) 


jLuminosity (pb~^) 


22.0 (60.0) 


116 (340) 


540 (1500) 


Channel: £-f}r 2j 


mg (GeV) 


300 


600 


900 


Signal events 


12.7 (35.3) 


7.60(21.1) 


6.48(18.0) 


Background events 


14.1(39.3) 


4.28(11.9) 


2.77(7.70) 


/Luminosity (pb^^) 


6.80(18.9) 


37.1(103) 


205 (570) 


Channel: i+fr 2j 


mg (GeV) 


300 


600 


900 


Signal events 


17.8(49.4) 


12.8(35.5) 


8.78(24.4) 


Background events 


29.7(82.5) 


14.1(39.1) 


6.08(16.9) 


jLuminosity (pb~^) 


11.2(31.1) 


75.6(210) 


283(785) 



discovery of heavy quarks corresponding to three different channels, namely, 2j, l~ fjr 2j 
and l^fjT 2j. While the number of signal and background events are model- independent 
quantities, the luminosity required is model-dependent. The luminosities shown are for a 
U,D and U type quark for the 2j, i~pr 2j and i~^f)r 2j channels respectively. For 
each of these quarks, the model-dependent parameter, k = 1 while the decay of the heavy 
quark in the respective gauge boson mode is 100%. The choice of these parameters is for 
illustrative purposes only and corresponds to the degenerate bidoublet scenario. 

We present the model-independent results of our analysis for the LHC in Fig. M for the 
different channels. There are three parameters - S"^, mg and the luminosity - in our analysis 
and we present our results in the planes defined by two of these parameters. In Fig. M^a) 
and Fig. MyC) we present the sensitivity plot for 3a evidence and 5a discovery in the plane 
of the parameters S^~, S'^'^ and and heavy quark mass tuq for i/s = 7 TeV with 1 fb~^ 
of data and for ^/s = 14 TeV with 100 fb~^ of data respectively. In Fig. Mjo) and Fig. M^d) 
we present the luminosity needed for 3a evidence and 5a discovery vs 5*^ at the LHC with 
y/s = 7 TeV and -y/i = 14 TeV respectively. In Figs. M (b) and (d) the quark masses change 
along the different curves. For reference we show with stars three different values of quark 
masses on each curve. 

The optimal channel for discovery in an early LHC run is rriQ ^ 1 TeV for 5*^ ~ 0.41pb 
for ^/s = 7 TeV and 1 fb~^ integrated luminosity. This is comparable to the final Tevatron 
reach for single production of these quarks. Of course, the reach increases dramatically for 
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TABLE VII: Same as Table EJ but for 14 TeV. 



Channel : 2j 


mq (GeV) 


900 


1800 


2400 


Signal events 


10.1 (28.3) 


8.60 (24.1) 


13.1 (37.0) 


Background events 


8.20 (23.0) 


5.80 (16.3) 


15.1 (42.8) 


jLuminosity (fb~^) 


0.15 (0.42) 


2.64 (7.40) 


35.0 (99.0) 


Channel: I- pr 2j 


mq (GeV) 


900 


1800 


3000 


Signal events 


9.29(25.8) 


4.72(13.1) 


4.18(11.6) 


Background events 


6.70(18.6) 


1.22 (3.40) 


0.720(2.00) 


jLuminosity (fb^^) 


0.0450(0.125) 


0.360(1.00) 


6.12(17.0) 


Channel: l+pr 2j 


mq (GeV) 


900 


1800 


3000 


Signal events 


15.7(43.7) 


13.1(36.3) 


6.48(18.0) 


Background events 


23.6(65.6) 


14.4(40.0) 


2.41(6.70) 


jLuminosity (fb^^) 


0.0864(0.240) 


1.15(3.20) 


11.5(32.0) 



a longer LHC run. We have estimated that with 100 fb"^ at 14 TeV and S"^ = 2.6 fb, the 
LHC can discover a new Q quark with a mass up to mq ~ 3.7 TeV. 

C. Constraints on Specific Models 

We interpret the model-independent results presented in Fig. [9] in terms of some specific 
models next. We present the sensitivity of the model dependent parameter kqQ as a function 
of the mass of the heavy quark (mg) for some popular scenarios. We present the 3o" and 
5cr constraints on kqQ for ^/i = 7 TeV with 1 fb~^ of data and a/s = 14 TeV with 100 fb~^ 
of data in Fig. [10] and Fig. [11] respectively. We also show the model dependent constraints 
from electroweak precision observables, which correspond to 

KqQ > 0.75, 0.066, 0.041, 0.037, (32) 

for the degenerate bidoublet, doublet, D singlet and U singlet models, respectively. We 
also show in Fig. [10] the current constraints from the recent Tevatron searches. Several 
conclusions can be obtained from our results. First, as expected, only for the degenerate 
bidoublet model can the Tevatron and the early LHC run probe regions of parameter space 
which are not excluded by electroweak precision tests. Second, the early LHC run reach is 
comparable and in general better than the Tevatron reach. In some cases, like the neutral 
channel in the U singlet model, the Tevatron is unable to constrain the model (the probed 
cross sections cannot be generated due to the bound from the unitarity of the mixing) 
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FIG. 9: (a) Top left: Sensitivity plot for Scr evidence and 5a discovery in the plane of the 
parameters S^~,S^^ and and heavy quark mass uiq with 1 fb~^ of data at ^/s = 7 TeV; 
(b) top right: luminosity needed for 3cj evidence and 5a discovery vs the parameters S^~,S^~^ 
and at the LHC with = 7 TeV. The three points on each curve represented by black (blue 
or red) stars correspond to the mass for each S^{S^ or S^^); (c) bottom left: same as (a) but 
for = 14 TeV and 100 fb^^ of data; (d) bottom right: same as (b) but for ^/s = 14 TeV. The 
vertical solid green and vertical solid dashed green bounds are the 95% C.L exclusion limits from 
Tevatron for the heavy quarks produced in pairs via QCD interactions and decaying via CC and 
NC decay modes respectively. 



whereas they are all probed by the LHC. Finally, the reach of the LHC run at y/s = 14 TeV 
and 100 fb~^ integrated luminosity is sensibly better, probing regions of parameter space 
not excluded by electroweak precision tests for all four models (for large enough masses). 
In particular, masses in the 3 — 4 TeV region can be probed in all four models provided we 
have kgQ > 0.4. Recall that for these large masses, kgQ can be directly related to a Yukawa 
coupling in particular models. Next we consider a few example models in detail. 
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1. D type singlet 

The first model contains an extra vector-like electroweak quark singlet with hypercharge 
— 1/3. We assume the singlet mixes only with the first generation (doublet) in the basis in 
which all SM flavour occurs in the up quark sector. The model can be fully parameterized 
in terms of the physical down quark mass, md-, the physical heavy down quark mass, mr,, 
and the coupling between D and m, k^^, = VudSi^ with V the CKM matrix and s^, is used to 
denote the sine of an angle (and is therefore smaller than one). The decoupling limit occurs 
for m^i — > oo. Si — i- with SLTfiD constant. In that limit all effects go to zero. The non-zero 
couplings apart from k^^,, are 

i^dD = slCl, (33) 

which is therefore smaller than 1/2, and Yukawa couplings given in the notation of Ap- 
pendix [B] 

YdD = slCl , Yod = slCl — . (34) 

V V 

These Yukawa couplings are needed to compute the corresponding branching fractions. 



2. U type singlet 

A very similar model consists of an extra singlet U with hypercharge 2/3 and we assume 
that it mixes only with the flrst generation (doublet) in the basis in which all SM flavour 
occurs in the down quark sector. The model can be fully parameterized in terms of the 
physical up quark mass, niu, the physical heavy up quark mass, mu, and the coupling 
between U and d, i^^^j = VudSi, with V the CKM matrix and sl is used to denote the sine 
of an angle (and is therefore smaller than one). The decoupling limit occurs for mu — )■ oo, 

— i- with slttiu constant. In that limit all effects go to zero. The non-zero couplings 
apart from k^^, are 

i^uu = SlCl, (35) 

which is therefore smaller than 1/2, and Yukawa couplings given in the notation of Ap- 
pendix [B] 

Yuu = SlCl , = slCl — • (36) 

V V 

These Yukawa couplings are needed to compute the corresponding branching fractions. 



3. Doublet with Hypercharge 1/6 

We now consider one hypercharge 1/6 doublet that only mixes with un in the basis in 
which all SM flavor occurs in the down sector. Everything can be parameterized in terms of 
the up quark mass, m„, the D quark mass, m/j, and k^d = —Sr- The other non-vanishing 
couplings are 

i^uu = -srCr, Hdu = -i^uD = Sl, (37) 
and a Yukawa coupling between u and U that we use to compute the branching fractions 

Yuu = srCr — , Yuu = srCr , (38) 

V V 
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where the mass of the heavy U quark is 



mu = ^ ^ ^ (39) 
The two LH couphngs are irrelevant in practice, since we have the relation 

sl = ~ 0, (40) 

for above current experimental limits. 



4- Degenerate Bidouhlet with Hypercharge 1/6 and 7/6 

This model discussed in detail in Ref. 0] and in Sec. [TTl consists of two new vector-hke 
quark electroweak doublets, with hypercharges 1/6 and 7/6 respectively, that mix only with 
the up quark, in the basis of diagonal Yukawa couplings for the charge 2/3 SM quarks. The 
two doublets are degenerate, i.e. they have the same mass and Yukawa couplings to u before 
EWSB, and the particle content of these two doublets is two charge 2/3, one charge —1/3 
and one charge 5/3 quark. The three quarks of types U, D and X in our notation, have the 
following values of k 

KuU = V2KuD = V2KuX = Sr, (41) 

where spt is the sine of the corresponding mixing angle which depends on the particular 
values of the model parameters. In this model, the generic bound from the unitarity of the 
mixing is saturated for Kuu < 1 and is a bit more stringent for k^^) — f^ux 

<1/V2. Also 

the constraints from the oblique parameters impose a mild bound of KuU ^ 0.75 (assuming 
AS < 0.2) in the model under consideration. 

In this model each of the quark decays with 100% branching ratio to the respective gauge 
boson, namely Br(f/ ^ uZ) = Bt{D uW) = Br(X ^ uW) = 100%. 



D. Systematic Effects 

The analysis reported here gives only an estimate of the LHC reach. Although better 
results could be obtained if the signal selection was based on a Likelihood or Discriminant 
Variable technique, it must be noted that a number of systematic effects, which are difficult 
to estimate by simulation, have been neglected. In particular, 

• No k-factors have been taken into account for the signals and backgrounds. 

• Since the charge identification plays a significant role in the CC analysis and Delphes 
does not simulate charge misidentification, the likelihood of charge mismeasurement 
was studied at the generator level. As an example, pr and rj dependent parameteriza- 
tion of the ATLAS tracker response was taken from j45t]. We find the misidentification 
fraction to be very small, only going up to 1.2% (0.9%) for the electrons (muons) from 
the heaviest quark (3TeV) studied in this paper. 

• Multijet events, which are abundant at the LHC, can fake the CC decay signal if jets 
are misreconstructed as leptons and if jet transverse energies are poorly measured, 
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FIG. 10: (a) Top left: sensitivity plot for 3a evidence and 5a discovery in the plane of the model 
parameters KuU and heavy quark mass niQ at the LHC with 1 fb^^ of data and ^/s = 7 TeV for 
the model with degenerate bidoublets; (b) top right: same as (a) but for kuD for the model with 
one doublet; (c) bottom left: same as (a) but for kuD for the model with D type singlet (d) bottom 
right: same as (a) but for k^u for the model with U type singlet. The curves with data points 
show the 95% C.L. constraints on kgQ for each model coming from the direct searches by CDF and 
DO experiments in the channel where heavy quarks are produced singly via EW interactions. The 
region above the blue solid line denoted "EW Constraints" is excluded by the EW precision tests 
described in section IIIBI The bounds on the other couplings in each model can be obtained from 
the relations between couplings listed in Sec. IVI Ci 



leading to a presence of missing energy in the events. These effects are very difficult 
to estimate by fast Monte Carlo simulation since they depend strongly on the criteria 
used for lepton identification and jet measurements. Nevertheless, we expect that with 
optimized experimental methods this source of background will be small. Indeed, for 
the DO and CDF analyses 0, Si it is almost negligible. 
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FIG. 11: Same as Fig. lO] but for ^=14. TeV and 100 fb^i of data. 

• The selection cuts have been optimized for specific mass values. In practice, a sliding 
window search would have these cuts re-optimized at intermediate mass points. In 
some cases, because of experimental uncertainties and to have more robust results, 
the cuts could be loosened. For example, the p^"^ cut on the tag jet plus heavy quark 
candidate system, as given in Eq. (1271) . could be raised to 150 GeV for the 14 TeV 
analysis, at the cost of ~ 10% loss in significance. 

• At high luminosity, pileup effects can become very important. Pileup will affect the 
purity of the signal. The isolation efficiency and resolution of leptons will deteriorate 
and fake forward jets would be generated. 

VII. SUMMARY 

The Fermilab Tevatron has been leading the way for searches for new heavy particles, and 
the historical mission of the LHC at the energy and luminosity frontier will take the field of 
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high energy physics to a new era. Among many highly motivated new particles in the TeV 
region, are new vector-like quarks that can have sizable couplings to first generation quarks 
without conflicting with current experimental constraints. The coupling with valence quarks 
and unique kinematics make single production at hadron colliders the optimal discovery 
process, potentially superior to the pair production mechanism via QCD interactions. 

In this article, we first motivated the consideration for the new heavy quarks by presenting 
a few theoretical models. We then discussed the existing experimental bounds on the quark 
mass and their couplings to the SM fields. In order to make our study as widely appUcable 
as possible, we have performed a model independent-analysis of the discovery reach at 
the Large Hadron Collider for new vector-like quarks considering single production and 
subsequent decays via electroweak interactions. After an exploration of the unique features 
of the signal kinematics, we performed an exhaustive optimization for signal extraction over 
the SM backgrounds. With this optimized analysis, we determined the reach for the signal 
as a function of the heavy quark mass. The result is given in terms of model independent 
quantities that can be identified with production cross section times branching ratio of the 
heavy quark decay to the corresponding electroweak gauge boson and a SM quark. We then 
translated these model independent bounds into bounds on the parameters of four popular 
models, with two degenerate doublets of hypercharges 1/6 and 7/6; one single doublet with 
hypercharge 1/6; one down type singlet and one up type singlet. 

Wc considered two scenarios, one with a short LHC run with 1 fb^^ integrated luminosity 
at ^/s = 7 TeV and a longer run with 100 fb^^ at ^/s = 14 TeV. The reach of the former has 
been found to be comparable, although typically slightly better, to the final Tevatron reach 
with 10 fb~^ integrated luminosity. It is only in the case of the degenerate doublets the 
probed region by the Tevatron and the short LHC run are unconstrained by precision data. 
However, the reach dramatically improves both in terms of mass and couplings for the longer 
run, probing in all four models considerable regions of parameter space not constrained by 
electroweak precision observables . In particular, masses up to ~ 3.7 TeV can be probed 
with couplings of order one or equivalently, production cross sections times branching ratios 
as small as 2.6 fb. It is evident from these results that the discovery potential for new heavy 
quarks at the LHC is very encouraging. 
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Appendix A: Explicit calculation of S 



The one loop contribution to the S parameter from an arbitrary number of fermions is 
given by joj 



(Al) 



where the sum runs over all fermions, is the physical (diagonal) mass matrix, [/, 



are the couplings of left and right handed fermions to W^, respectively and Y^' the cor- 
responding couplings to B^, all in the physical basis. We have also defined the following 
functions 



x+iyi,y2) 



5{yt + yl) - 22y\yl ^ 2yiyl{y\ + y^) 



X- (1/1,1/2) 



log 

V\V2 



vl? 



log 7^ 



y± 



+ log 



(2/1 - yl 



2N3 



yt 



yl-2yiyl\og[^-\ 



(A2) 
(A3) 



The model with two exactly degenerate doublets that only mix with the up quark can be 
parameterized by three parameters which are usually taken to be the up Yukawa coupling, 
the Yukawa coupling between the up and the two doublets and the (common) mass of the 
two doublets. All these parameters are defined in the current eigenstate bases, i.e. before 
diagonalizing the mass matrix. It is however much more advantageous to define physical 
parameters which, in our case can be taken to be the (physical) mass of the up quark m„, 
the physical mass of the three degenerate heavy states (which coincides with M) and the 
coupling Kyxj = sr, which coincides with the sine of the RH rotation needed to diagonalize 
the mass matrix. Full details of this parameterization can be found in Ref. [l8|. The physical 
masses, written in terms of these parameters and in a basis with first the (three) charge 2/3 
quarks, then the two charge —1/3 quarks and finally the charge —5/3 one are. 



M = diag(m„, M, 



2^2 



Cr 



md,M, M), 



(A4) 



whereas the couplings relevant for the calculation of the S parameter are, in the same basis 
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(A6) 



where the LH rotation is given in terms of our input parameters by s/, = SRirnu/M (see 18|). 
Eq. flA6P shows that we indeed have Kuu = sr. We can now insert the exphcit form of the 
couphngs and masses in Eq. flAip and expand in the small ratio m^/M < 10~^. The resut 
is 



S = S^ 



SM 



64 - 94 + 54 + 3(2 - 44 + 34) log[l - s] 
3^4 



where 



SM 



3-log(m^/m^) 



(A7) 



(A8) 



is the SM contribution from the up and down quarks that we have to subtract to compute 
the new physics contribution to the S parameter. Thus, we see that, up to tiny corrections 
^('^u/-^^) ^ 10"^^, the new physics contribution to S in this model is independent of M 
for fixed values of Kuu- 



Appendix B: Heavy fermion decay widths 



In this appendix we collect the general expressions for the decay widths needed in the 
calculation of the branching ratios for arbitrary models. The decay width for a two body 
decay of a heavy fermion, Q with mass mg, into a light fermion, q with mass and a 
massive gauge boson, V with mass my, and couplings 



V.,, 



qi^igiVL + 9rVr)Q + h.c. 



reads 



T{Q ^ Vq) 



\Jra'^ — 2m^(mQ + niy) + (mg — niyY 



X \i\9L\' + \gRmml 



m 



2\2 



+ myimq + — 2m' 



-^{9l9*r + gl9R)mQmqml 



327r m 



V 



1 



ml 



m 



+ 2- 



m 



V 



Q 



mV 



Similarly the decay into a quark and a scalar, with couplings, 

1 



V2 



H 



yqqQhQR + yqqqiQR + h.c. 
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(Bl) 



, for {mQ,mv > rrig), (B2) 



(B3) 



is given by 



r(g ^ Hq) = 



^JmjJ + {ml - - 2m|^(m| + m^) 



X 



647rmQ 



(B4) 



^Qq + ^gQ 

647r 



1-2- 



H 



+ 



m 



H 



, for {mQ,mH > mj. 
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